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Abstract 

Results of relativistic calculations of the direct knockout (DKO) 
mechanism for the photon induced removal of a proton from a target 
nucleus over a wide range of energies and nuclei are presented. Spec- 
troscopic factors used in the calculations are fixed from consistent 
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analyses of the quasifree electron scattering process (e, e'p). The re- 
sults indicate that within the uncertainties of the model, the knockout 
contributions are generally close to the experimental data for missing 
momenta below ~ 500 MeV/c. This is in disagreement with nonrcla- 
tivistic analyses which often find that the direct knockout contribution 
can be quite small compared to the data and that meson exchange 
corrections can be important. The present study suggests that meson 
exchange current contributions may not be as large when treated in a 
relativistic framework. We also point out some difficulties we encoun- 
tered in analyzing the data for a ^^C target at photon energies below 
80 MeV. 
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1 Introduction 



The reaction mechanism leading to the knockout of a single proton by a real 
photon has been the subject of some debate recently. Some nonrelativistic 
analyses JI], 0, 0, |[ suggest that the direct knockout (DKO) contribution 
may be very small compared to the data. Accordingly it was concluded 
that meson exchange current (MEC) contributions must be the main mech- 
anism responsible for the observed cross sections. Similar conclusions were 
reported earlier for the nonrelativistic analyses carried out by Miller et al. 
[§] for ground state transitions for the reaction ^^C'(7,p)^^N at 60 and 72 
MeV, and by Ireland et al. for (7,p) reactions on several nuclei for 
near 60 MeV (the same nuclei involved in the discussion reported in Ref. 
[H]). The above conclusions do not seem consistent with the findings by 
Ryckebusch et al. 0. These authors find MEC effects to be relatively small 
for ground state transitions. The above statements illustrate the existing 
difficulty of arriving at a consensus within the nonrelativistic framework as 
to the extent of contributions from processes beyond simple direct knock- 
out to ground state transitions in (7,^) reactions. Although the differing 
views stated above appear to be somewhat dependent on the nuclear models 
used in the nonrelativistic calculations, they are, however, symptomatic of 
our incomplete understanding of the nature of the reaction mechanism for 
photonuclear reactions. 

These results are quite different from those of a recent relativistic analysis 
of Johansson et al. who find that, for an incident photon energy of 60 
MeV, the DKO contribution accounts for most of the observed data, with no 
indication of any systematic sharp deviation from the data at this energy. 

In this paper we extend the analysis reported in Ref. [Q, for a photon 
energy of 60 MeV, to a much wider range of data. We consider several data 
sets for the (7,p) reaction [and some data on the inverse reaction (p, 7)] on 
a number of target nuclei and covering a range of photon energies extending 
well into the A-resonance region. The spectroscopic factors and wave func- 
tions used in the calculations are fixed at the values obtained from a parallel 
analysis of the (e, e'p) reaction on the same target nuclei. The objective of 
this study is to use this type of constrained analysis to gain some insight into 
the role of the DKO mechanism and to see if a consistent description of the 
available data is possible. 

Section || outlines the relativistic calculations for the direct knockout 
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contribution to the (7,^) reaction. Results of the calculations and details of 
the comparisons with data are given in section ^. Our conclusions are given 
in section ^. 



2 Relativistic Calculations 



The differential cross section due to the direct knockout contribution to the 
(7, p) reaction has been given previously P, |^ but we provide it here again for 
ease of reference. The relativistic expression for the differential cross section 
leading to a specific final state of the residual nucleus can be written as 



da 

dQr, 



a Mc^ 
Air he 



Pp 



C C 1 Sj^Jj{Jb) ^ L/3;yM/ 
Vrel R 2 Jb + 1 ^^^^ ^ 



(1) 



where Mb and /i are the spin projections of the bound and continuum pro- 
tons. We denote the four-momentum of the final proton pp and the four- 
momentum of the incident photon as q. The four-vector is the photon 
polarization vector with two polarization states r, and summation is implied 
over repeated greek indices. The recoil factor R is given in any frame by |]10| 



R = 1 



E 



R 



2Pp ■ Pr- 



(2) 



The four-momentum of the recoil nucleus is denoted by pr. The function 



N^^'' is 



J d^x {pp, x) r/3^js,Ms (x) exp (iq ■ x) , 



(3) 



where the wave functions of the continuum and bound nucleons, denoted 
and "^jg^Ms respectively, are solutions of the Dirac equation containing 
appropriate potentials The 4x4 matrix F^, operating on the nucleon 
spinors, is given by 



/3 = 7o 



7/3 + 



2M 



P u 



(4) 



The ingredients of the model are basically the same as used in an analysis 
of light to medium weight nuclei at 60 MeV [Bl: the bound state protons are 
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described by solutions of a Dirac equation containing the relativistic Hartree 



potentials of Blunden and Iqbal [jTT| , while the final state continuum proton 
is described by solutions of a Dirac equation containing complex phenomeno- 
logical optical potentials obtained from fits to proton elastic scattering data 



1^. Given these potentials, the only parameters left to determine are the 
spectroscopic factors. For the light nuclei: ^°B, ^^C and ^^O, we have ob- 
tained the spectroscopic factors by fitting the results of our (e, e'p) model 



[§, |T3[ to available data. The ^°^Pb data are not suitable for analysis using 



this model because of the lack of Coulomb distortions for the incident and 
final electrons. The spectroscopic factors used in this case are those of Udias 
et al. |14| who have performed a relativistic analysis of the ^osp^ data. 



In the following we show the results of our calculations compared to the 
experimentally determined cross sections for several nuclei covering a wide 
energy range. 



3 Results and Discussion 

We have performed calculations for the (7, p) reaction on several target nuclei, 
over a wide range of energies. These are compared to existing data in order 
to assess the extent to which the direct knockout mechanism contributes to 
the observed cross sections. The ingredients of the calculations have all been 
determined elsewhere and since there are no adjustments made, the results 
can be considered as predictions of the model. 

In the graphs to be discussed below there are curves corresponding to 
several different calculations. The description of the calculations represented 
by each of these curves is as follows: 

1. dashed curve — energy- (E-) dependent parameterization of the Dirac 



optical potentials specific to a single nucleus while the bound state 



wave function is obtained through a Dirac-Hartree calculation P]] ; 

dotted curve — E-dependent parameterization of the Dirac optical po- 
tentials specific to a single nucleus and the binding potential has a 
Woods- Saxon form; 

solid curve — energy- and mass- [(E+A)-] dependent parameterization 
of the Dirac optical potentials and the same Dirac-Hartree bound state 
wave function as in curve 1 above; 
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4. dot-dashed curve — curve 3 divided by a factor of 2.0 to bring the 
model calculations close to the data. 

All the figures shown below use this designation of curves. The first three are 
simply for calculations using a variety of existing potential models in order 
to provide some feeling for the sensitivity of the results to variations in these 
ingredients. The dot-dashed curve is only relevant to graphs shown for the 
^^C target. 



3.1 12 c Target 

A considerable amount of data are available for this target. We have made 
comparisons of our relativistic DKO model calculations with these data, con- 
centrating mainly on ground state transitions. These comparisons are shown 
in Figs. 1-3. Data for four of the energies shown in Fig. 1: = 49.0, 58.4, 



67.8 and 78.5 MeV were reported by Springham et al. ||T5| and are obtained 
using the tagged photon facility at Mainz. The absolute magnitude of their 
cross sections was obtained by normalizing the data at each energy to data 
taken by Mathews et al. 0] for the (7,^0+1) reaction, data which include 
both the ground and first excited state of the residual ^^B nucleus. The data 



of Aschenauer et al. [|T7|, at photon energies of = 45.0 and 54.0 MeV 
were obtained at the MAX-Lab at the University of Lund. These data were 
normalized completely independently of any previous experiment, and found 
to be consistent with existing data within systematic errors. The data at E, 



1 



73.5 MeV from Rauf [181 were also obtained at the MAX-Lab. These data 



were normalized to previous measurements including those of Mathews et al. 



le 



The most obvious feature apparent in Fig. 1 is that the calculations tend 
to form a narrow band lying above the data. Note that there is not much 
sensitivity to reasonable variation of the ingredients of the model. The light 
dot-dashed curve in all the figures shows the solid curve divided by a factor 
of 2, and this brings the curve close to the data in all cases. The calculated 
curves have the correct shapes and the variation of magnitude with incident 
photon energy also seems to be correct, but the curves lie consistently above 
the data by a factor of 2. 

Figure 2 shows the differential cross section as a function of photon en- 
ergy for four different proton angles. The experimental data are taken from 
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Ruijter et al. [|T9|; the experiments were also performed at the MAX-Lab at 
the University of Lund. Absolute normalization of these data was obtained 
independent of any other experiment and the results were found to be con- 
sistent with a large amount of other data. Again the calculated curves lie 
above the data by close to a factor of 2. We see that for photon energies 
above 40 MeV the energy and angular dependence are quite well reproduced 
by our model, but results lie above the data by a factor of 2. 

In an attempt to compare to other data as well, we consider two exper- 
iments in which the first excited state of ^^B at 2.12 MeV (1/2^) could not 
be resolved from the ground state. The experiment of Mori et al. ^ was 
performed at the Laboratory of Nuclear Science at Tohoku University. The 
experiment of Harty et al. was performed at Mainz. The differential 
cross sections of both experiments were normalized without reference to any 
other experimental results and found to be consistent with other data. In 
order to compare to some data which do not contain the first excited state, 
ground state data from other sources were multiplied by a factor of 1.27 by 
both groups. This factor is assumed to account for population of the first 
excited state in ^^B being ~27% as probable as population of the ground 
state over the range of energies considered in both experiments. We adopt 
this factor in what follows. 

In the left-hand column of Fig. 3 we show the measured energy distri- 
butions of Mori et al. B as the circular data points. The triangular points 



are taken from the data of Harty et al. ||2^. The angles at which these data 
were taken do not coincide with the angles of the experiment by Mori et al. 
For this reason the triangles shown at 66.0° represent cross sections which 
have been averaged for proton angles of 63.3° and 68.4°, while on the graph 
labeled 91.1° we show data averaged for 88.5° and 93.5°. Note that the 88 
MeV data point of Harty et al. lies almost on top of the 87.8 MeV data 
point of Mori et al. showing the consistency between the two data sets. The 
curves are calculated assuming that the recoil nucleus is left in its ground 
state, and then multiplied by 1.27. The curves again lie above the data by 
a factor of 2 at low photon energies, but at the higher energies of the Mainz 
experiment the calculations seem to move closer to the data. 

The right-hand column of Fig. 3 shows the angular distributions obtained 
by Harty et al. compared to calculations as discussed in the previous 
paragraph. As the energy of the incident photons increases, the calculations 
seem to move from lying above the data by about a factor of 2, to falling 
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within the error bars for the data at the highest energy. Of course, because 
of the large error bars at larger proton angles the trend is not definitive, but 
it is suggestive. 

The picture that is emerging for the status of the comparisons for the 
^^C case can be further clarified by looking at the data available for {p, 7) 
reactions on this nucleus as well as those leading to its formation as a resid- 



ual nucleus. The data are those obtained recently by Bright et al. |£T1 at 
Uppsala. Figure 4 shows comparisons to the data at proton energies of 98 
and 176 MeV for ground state transitions to ^^C and ^^N residual nuclei. 
The former reaction is the inverse of the (7,p) reactions discussed above. 
The comparisons for the ^^C residual nucleus are shown on the right hand 
side of the figure. Using the same wave functions and spectroscopic factors 
as in Fig. 1, we find that the calculations for Tp = 98 MeV lie slightly above 
the data at all angles except for the last point at 6p =140°. At 176 MeV 
the calculations are closer to the data except for the large angles. Thus the 
data for the inverse reaction confirm the behavior alluded to above; at lower 
energies the calculations seem to overestimate the cross sections. 

The comparisons on the left-hand side of Fig. 4 present a somewhat 
different picture. Using the maximum value for the spectroscopic factor, the 
relativistic calculations for radiative capture on ^^C are close to or below the 
data. With a more realistic value of the spectroscopic factor the calculations 
will be further reduced in magnitude. This situation is in clear contrast to the 
cases discussed above. It should be noted, however, that with a spectroscopic 
factor in the range 0.5-1.0 (the maximum possible value is 1.0 in this case), 
the contributions of the knockout mechanism to the reaction are substantial 
in the region of lower missing momenta. 

It is worthwhile to point out here that there are also unresolved difficulties 
for the ^^C target in the (e, e'p) reaction, in addition to the difficulties dis- 



cussed in the current work. The data from NIKHEF |22[ are for kinematics 
with a fixed final proton kinetic energy of Tp = 70 MeV and nonrelativistic 
calculations shown in that paper cannot reproduce the shape of the distri- 
bution in missing momentum. In particular when the calculations are scaled 
to fit the peak for positive missing momenta, the calculations fall below the 
data for negative missing momenta. Our relativistic calculations show ex- 
actly this behaviour and the spectroscopic factor that we have used in this 
work is obtained by matching to the positive missing momentum peak. A 
proposed solution to this problem was to adjust the ratio of transverse to 
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longitudinal response functions, and when this ratio was adjusted to ^1.3 
the shape of the missing momentum distribution was reproduced. This 
problem was considered further by van der Steenhoven who found no 
justification for this enhancement factor. The newer data from Mainz on 
this nucleus, reported by Blomqvist et al. p4| have a higher final proton ki- 
netic energy of Tp ^ 90 MeV. In this case, with the increase in normalization 



of the data by a factor of 1.19 both nonrelativistic and our relativis- 
tic calculations can describe the shape of the measured missing momentum 
distribution using the spectroscopic factor as obtained from matching to the 
positive missing momentum peak of the NIKHEF data. This behaviour is 
consistent with the current results for the (7,p) reaction on this nucleus: at 
low final proton energies model calculations differ from the data and as the 
final proton energy increases the calculations move closer to experimental 
results. 

The data for proton knockout from the lp3/2 orbital in ^^C seem to in- 
dicate that at low energies a simple shell model description is not adequate 
to explain the data. A more complete description might possibly involve 
inclusion of the deformed nature of the ground state wave function through 
a configuration mixing picture. 



3.2 ^"B Target 

Data for both the {e,e'p) and (7,p) reactions has been obtained for a 
target by de Bever p6|] at two different energies. The data for knockout of a 
1^3/2 proton leading to the ground state of ^Be are shown in Fig. 5. Spectro- 
scopic factors were obtained by scaling the model calculations to the (e, e'p) 
momentum distribution (or reduced cross section) data for the Tp = 70 MeV 
case. The other (e, e'jo) and (7,p) curves were then calculated without any 
adjustment of the parameters. It should be noted that the optical potentials 
used here are parameterized using proton elastic scattering data on targets 
from ^^C to ^°^Pb. As a result of a lack of proton elastic scattering data 
on ^°B, we simply extrapolate the (E+A)-dependent potentials for use with 
a target hghter than ^^C. The {j,p) calculations for this nucleus are quite 
sensitive to changes in the potentials used to generate the nuclear wave func- 
tions. In spite of this, and the fact that ^"B is not a closed shell nucleus, the 
DKO clearly produces results in the neighbourhood of the data. 
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3.3 16 O Target 

Figure 6 shows the differential cross section as a function of proton angle 
for knockout of a valence proton from an ^^O target, leading to the ground 
state of ^^N. The photon energy range is the largest available, with eight 
energies in the range 60 MeV < < 361 MeV. The data come from three 
sources: Miller et al. |^ provide data points at energies of 60 an 72 MeV, 



while data shown at 60, 80 and 100 MeV are from Findlay and Owens 
The high energy data for photons in the range 196 MeV < < 361 MeV are 



from Adams et al. At low energies the calculated curves are generally 



close to the data, reproducing the magnitude and shapes quite well. For the 
higher energy data of Adams et al. the calculations tend to be close to the 
data points at small angles while falhng below the data as the proton angle 
increases. This is the behavior one expects if meson exchange processes are 
going to become important as the missing momentum increases. 

In order to remove some of the kinematic dependence from these curves 
we have calculated a reduced cross section by dividing the differential cross 
section of Eq. (|ip by a kinematic factor [^, [I|]: 



Ep 



E^ M2 



\m'{ep) + ^KlE'^ 



(5) 



Figure 7 shows the reduced cross section as a function of missing mo- 
mentum for all the experimental data shown in Fig. 6, as well as additional 
data provided by Leitch et al. |]5D[ . The curves are generated using the same 



ingredients as the solid curves of Fig. 6 but restricted to the kinematic range 
covered by the data. An interesting observation here is that the model re- 
sults are close to the data for missing momentum less than about 500 MeV/c. 
The vertical dotted line indicates the momentum of a free proton with ki- 
netic energy equal to the charged pion mass. The calculations start to fall 
below the data in this kinematic region, which seems to be a good indication 
that we are seeing the need for inclusion of pion exchange diagrams to the 
reaction mechanism, and provides some idea of where these diagrams become 
important. 
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3.4 ^"^Pb Target 

Figure 8 shows results for proton removal from ^^'^Pb, leading to two doublets 
and one resolved state in ^°'^T1, for two relatively low photon energies: 45 
MeV and 54 MeV. The data are from Bobeldijk et al. These authors 
performed a nonrelativistic distorted-wave impulse approximation (DWIA) 
analysis of the data and found that the DKO contribution tends to lie up to 
a factor of 10 below the data. Revised recent analyses |l^, ISTJ indicate that 
this factor may have been unrealistic. Our present analysis, on the other 
hand, shows that the relativistic calculations do come close to predicting the 
correct magnitudes of the observed cross sections. Within the parameter un- 
certainties, it is evident that the DKO mechanism is the leading contributor 
to the reaction at these energies. 

4 Conclusions 

In this paper we have presented relativistic calculations for the (7,^) reaction 
and its inverse for a number of target nuclei. The results for the light targets 
cover a wide energy range, while the results for the lead target are at low 
energy but for a variety of final nuclear states. The analysis was done in a 
consistent manner with no free parameters. In all cases but one, ^^C(p, 7), the 
spectroscopic factor is obtained from a parallel analysis of the corresponding 
(e, e'p) data. 

In cases of transitions with simple nuclear structure, relativistic calcula- 
tions indicate that the DKO mechanism is the main contributor to the cross 
section for lower missing momenta. For larger missing momenta one finds 
clear deviations indicating an increased role for higher order processes such 
as meson exchange and A-isobar contributions. 

Nonrelativistic analyses often indicate that the contributions from the 
DKO mechanism are small and that meson exchange effects are sometimes 
dominant even at lower energies. In contrast, the present relativistic analysis 
suggests substantial contributions from the DKO mechanism to the cross 
sections over a wide range of energies. The analysis also points out that 
meson exchange effects are required at higher missing momenta. 

In the course of this analysis we have found that in the case of the ^^C 
target for photon energies below 80 MeV, the relativistic calculations appear 
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to overestimate the cross section data by close to a factor of 2. This situation 
is puzzhng and may indicate either some comphcations due to the structure 
of the ^^C nucleus itself or to some subtleties in the combined analysis of 
(e, e'p) and (7,p) for this target. It must also be noted that these difficulties 
do not occur for the spherical nuclei ^^O or ^"^Pb, which are also considered 
in the present work. It is our feeling that the differences between theory and 
experiment at the lower proton energies for the ^^C target reflect the need for 
a proper description of the structure of this nucleus to include the intrinsic 
ground state deformation. The consistent approach based on a combined 
analysis of these two reactions |3^, |1|, || leads, in our view, to the conclusion 
that the ^^C ground state cannot be adequately described by simple single 
particle conflgurations. 

In the case of transitions with simple structure (mainly single particle) 
our calculations indicate that meson exchange effects will not be important 
until one reaches missing momentum near 500 MeV/c. With the effort to 
push (e, e'p) reactions towards this region of missing momentum it would 
be interesting to see how important MEC effects will turn out to be in the 
relativistic model. Van der Sluys et al. have considered this question in 
a nonrelativistic random phase approximation (RPA) framework and found 
large contributions from MEC's for larger missing momenta. 

One point of interest is that the reactions discussed, (e, e'p) and (7,p), 
show different sensitivities to the description of the bound state. This is 
probably due to the different range of missing momenta sampled by the two 
reactions. The (e, e'p) reaction has been primarily concerned with low missing 
momenta where the bound wave function is constrained by properties such as 
binding energy and rms radius. The bound state wave functions that we use 
show little difference in momentum space for small momenta, and so it is not 
surprising that the (e, e'p) results are very similar in this region. Differences 
between the bound state wave functions do arise, however, for larger missing 
momenta in (e, e'p) and for the inherently large missing momentum reaction 
(7,p). This is not a surprise because this is where the nuclear wave function 
is poorly constrained and the region where we see differences between these 
bound state wave functions in momentum space. 

A common criticism of the distorted-wave Born approximation (DWBA) 
approach, both relativistic and nonrelativistic, is the lack of orthogonality 
of the bound and continuum wave functions. It is argued that this lack of 
orthogonality could invoke spurious contributions to the cross sections. The 
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distorted continuum wave function is an approximation to the many-body 
wave function of the nuclear system with appropriate boundary conditions. 
This approximation derives its support from the fact that the wave function 
is constrained by proton-nucleus elastic scattering data. Nonrelativistic RPA 
calculations do not suffer from this lack of orthogonality, but the wave func- 
tions are not able to account for the elastic scattering data. A simple method 
for restoring orthogonality has been suggested by Boffi et al. [0] and Ciofi 
Degli Atti et al. These authors find the orthogonality effects to be 

relevant mainly at large angles. It is likely that this feature will carry over 
into the relativistic calculations and hence would not substantially change 
the main characteristics of the present calculations. 

The present results pose certain challenges for the relativistic approach. If 
the DKO contributions are large then the data would suggest that the MEC 
effects are suppressed in the relativistic models, at least at the lower energies. 
Relativistic models must explain this suppression and in the meantime face 
the challenge of accounting for the observed relatively large photoneutron 
cross sections. 

Spin-dependent observables are likely to play an important role in clari- 
fying the reaction mechanisms. It should be noted here that the cross section 
angular distributions for (7,^) reactions do not have much structure in most 
cases. The differences between competing models are then mainly differences 
in magnitudes, and hence may be related to normalization uncertainties in 
the models. When we discuss spin-dependent observables these normaliza- 
tion uncertainties cancel out and hence a better test of the model is likely to 
result. 
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Figure Captions 



FIG. 1. Knockout of a lp3/2 proton from a ^^C target leading to the ^^B 
ground state. Angular distributions for seven different photon energies rang- 



ing from 45 to 78.5 MeV. Hartree bound state wave functions are used [11 



and the proton optical potentials are from Ref . ||T2[ . The data are from Refs. 



1^, [jT^ and [|T8|. Curves as discussed in the text. 

FIG. 2. Knockout of a lp3/2 proton from a ^^C target leading to the ^^B 
ground state. Distributions in photon energy at four fixed proton angles: 6p 
= 30.0°, 60.0°, 90.0° and 120.0°. The data are from Ref. |l|. Curves as 
discussed in the text. 

FIG. 3. Knockout of a 1^3/2 proton from a ^^C target leading to the ^^B 
ground and first excited states. Left-hand column — distributions in photon 
energy at four fixed proton angles: 6p = 30.6°, 45.8°, 66.0° and 91.1°. Right- 
hand column — angular distributions for five photon energies The data are 
from Refs. and 0]. Curves as discussed in the text. 

FIG. 4. Differential cross section as a function of photon angle for the {p, 7) 
reaction on ^^C and ^^B leading to the ground state. Curves as discussed in 



the text. The data are from Bright et al. |21 



FIG. 5. Single proton removal from the ^°B target leading to the ground 
state in ^Be: upper figure - the (e, e'p) reaction, lower figure - the (7,^) 
reaction. Curves as discussed in the text. The data are from de Bever [^. 

FIG. 6. Differential cross section as a function of proton angle for the knock- 
out of a lpi/2 proton from a ^^O target leading to the ^^N ground state. 



Curves as discussed in the text. The dat from Refs. i, H P 



FIG. 7. Reduced cross section as a function of missing momentum for the 
knockout of a lpi/2 proton from a ^^O target leading to the ^^N ground state. 



Data as in Fig. 6 and from [Q. Curves as discussed in the text. 



FIG. 8. Differential cross section as a function of proton angle for the knock- 
out of protons from different levels in a ^^'^Pb target. Curves as discussed in 
the text. The data are from Bobeldijk et al. 
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